RADIO SCIENCE Journal of Research NBS/USNC-URSI 
Vol. 69D, No. 2, February 1965 

On the Point of Emergence of a Microwave Beam 
Entering a Linearly Graded Plasma 

A. L. Cullen 

Contribution From the Department of Electronic and Electrical Engineering, University of Sheffield, England 

(Received August 31, 1964) 

A beam of microwave radiation is assumed to enter a linearly graded plane-stratified 
plasma. If the electron density reaches the critical value cchotn/e 2 at a distance x c from the 
plasma boundary, ray theory predicts that the beam will emerge from the plasma at a 
distance 2x c sin 20 from its point of entry, being the angle of incidence. The validity of 
this result is studied using wave theory. 



1. Introduction 

The electron-density profile in a plasma has been 
studied by Kharadly [1963] using a microwave 
beam obliquely incident on the plasma boundary. 
By studying the dependence of the position at which 
the beam emerges again from the plasma as a func- 
tion of the angle of incidence, the electron-density 
profile can be calculated. The calculation assumes 
that ray theory is valid, and it is the purpose of the 
present paper to study the validity of ray theory 
for the special case of a linear variation of electron 
density with position. The ray theory result is 
very simple; the distance between the entering and 
emerging rays is 2x c sin 20, where x c is the depth in 
the plasma at which the electron density reaches the 
critical value at which the refractive index is zero, 
and 6 is the angle of incidence. 

2. Wave Theory of Reflection 

The propagation of electromagnetic waves in 
stratified plasmas has received a great deal of atten- 
tion in the literature, and the resulting body of 
knowledge has recently been summarized very con- 
veniently by Budden [1961], Wait [1962], and Ginz- 
burg [1964]. 

In particular, it is well known that electromag- 
netic wave propagation in a linearly graded plasma 
can be expressed in terms of Airy functions [e.g., 
Wait, 1962]. It is therefore possible to obtain an 
exact solution in closed form for the reflection 
coefficient for a plane wave in free space incident on 
such a plasma. When the electric field is perpen- 
dicular to the plane of incidence, the reflection 
coefficient, using Budden's notation [1961], is 
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I n the absence of collisions, k and f are real 
numbers, so 1/^1 = 1, and we can write R=e j *, where 
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If we consider a narrow angular spectrum of plane 
waves incident on the plasma, the reflected spectrum 
will be modified because of the dependence of \p on 
through (2) and (3). 

We wish to find the position at which the electric 
field of the reflected wave has a maximum value; 
let Ay be the distance (measured along the plasma 
boundary) between this position and the corre- 
sponding position for the incident wave. The 
principle of stationary phase can be invoked to show 
that 
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Using (2), we can write 



df__ 2/'(fl) 
dO l+j\B) 



(4) 



(5) 



The denominator of this expression can be written 
in terms of the Airy functions thus 

The calculation of the derivative f(6) is straight- 
forward, but is simplified considerably by using the 
result -4i // (fo) = fo^(fo) which follows from the 
differential equation which the Aiiy functions 
satisfy. Thus, we get 
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Combining (4), (5), (6), and (7), and rearranging 
the resulting expression with the help of (3), we 
finally get 
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It is convenient to use a parameter <p rather than 
fo, where 
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We can then define a correction factor F(<t>) by the 
following equation: 
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7^^) is therefore the factor by which the ray- 
theory result must be multiplied to get the wave- 
theory result. 



3. Asymptotic Formulas 

If (1) is simplified by using the first terms in the 
asymptotic expansions of Ai and Ai' , the reflection 
coefficient becomes 

R=j exp < — - j cos 3 OJcXc > (10a) 



so that 
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When (11) is substituted in (4), the ray-theory 
result is again obtained, thus verifying the well- 
known between ray theory and the W.K.B. approxi- 
mation. If the first two terms in each asymptotic 
series are retained, a first order correction to ray 
theory can be found, in terms of elementary functions 
and is as follows : 
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where <£>=<p— 7r/4. There is, however, little to be 
gained from the use of (12) rather than (8) as the 
Airy functions are quite adequately tabulated. 
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Figure 1. Ray theory correction factor. 

4. Discussion of Results 

In figure 1 the function F(<p) is plotted as a function 
of <p. The full line is the exact expression in (8); 
the broken line is obtained from (12) and shows how 
the asymptotic result fails for <£<37r/2. The Airy 
functions Ai and Ai f are tabulated [British Associa- 
tion Mathematical Tables, 1946]. 

It is clear that if <£>0.657r, ray theory and wave 
theory agree to within 5 percent. For example, 
if 0=60°, we must have # C >3.9A for 5 percent accu- 
racy. This condition is quite easily met in 
Kharadly's experiments, but it is important to bear 
this limitation of ray theory constantly in mind in 
applying the oblique incidence plasma diagnostic 
technique in the millimeter-wave region. A similar 
correction factor occurs in correcting ray theory in 
the vicinity of caustics [e.g., Wait, 1962, p. 338]. 
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